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PLASMA GENERATION USING A 
LARGE V-GROOVE CATHODE DISCHARGE 


By Paul R. Caron 
Electronics Research Center 


SUMMARY 

The results of an investigation of the plasma parameters 
generated by a large V-groove cathode discharge are presented. 

The cathode consists of a 12-inch aluminum plate into which 
V-grooves of 1/8-inch spacing and 60~deg angle are milled* With 
this cathode, a negative glow discharge is produced which yields 
a stable, steady-state plasma. Electron densities up to 2x10^2 
cm - ^ have been obtained in the 12- inch diameter chamber. The 
electron temperature is about 1200°K and is independent of dis- 
charge current and gas pressure over the range 0.3 to 1.0 torr. 

The discharge operates in much the same fashion as a con- 
ventional brush cathode discharge but with lower operating voltage 
for equal current densities. The lower operating voltage, which 
reduces the arcing problems at the cathode-glass junction, and 
the relative ease in fabricating the V-groove cathode make it an 
attractive alternative to the brush cathode for many applications. 
It is shown that the plasma is beam-generated and the dominant 
electron-loss mechanism in the steady-state is electron-ion 
recombination. However, although most electrons are lost by 
recombination, diffusion losses still play a strong role in 
determining the electron density profile. 


INTRODUCTION 

For many laboratory investigations, it is desired to con- 
struct a physically large plasma of sizable electron density. In 
particular, these characteristics are often desirable for study- 
ing microwave propagation in plasmas where the plasma should be 
many free-space wavelengths thick and where the propagation 
characteristics for wave frequencies near the plasma frequency 
are desired. The V-groove cathode discharge discribed in this 
technical note is well suited to such investigations. It will be 
shown that the plasma thickness and the electron density can be 
varied over a wide range. The plasma temperature is constant and 
the electron-neutral collision frequency can be varied by about a 
factor of 3; however, this cannot be done independently of the 
other quantities. Also, the electron-ion collision frequency is 
at least as large as the electron-neutral collision frequency. 



PHENOMENOLOGICAL DESCRIPTION OF THE DISCHARGE 


The discharge operates basically as a cold-cathode , abnormal, 
negative glow discharge in helium. The discharge configuration 
is shown in Figure 1. The tube diameter and, hence, the diameter 
of the cathode, Is approximately 12 inches* The anode- to-cathode 
distance can be varied at will but, for reasons which will be 
discussed later, the results of this report will be confined to 
the case wherein this distance is 6 inches* The cathode consists 
of an aluminum plate with milled V-grooves of 60-deg angle and 
1/8-inch spacing. With these grooves, the cathode current caused 
by secondary emission processes (i.e., ion-bombardment and photo- 
electric emission) remains uniform over the cathode surface In 
the abnormal glow region* That is, the tendency for a localized 
discharge or arc to form, which appears in large , flat cathodes , 
is reduced. Also, the grooves increase the effective cathode 
area and increase the electron yield by ion -bombardment due to 
non-normal incidence of ions on the cathode (ref. 1}* 

The chamber is pumped down to low pressure and backfilled 
with helium* To improve the gas purity, helium is allowed to 
flow continuously through the chamber at a slow rate and the pres- 
sure is adjusted by means of a leak valve in series with the 
helium supply. A dc voltage applied between the anode and cath- 
ode produces a negative glow discharge whose behavior is qualita- 
tively similar to that of a brush cathode discharge (ref. 2) . 

That is, the plasma is beam-generated and the electron density 
is sufficiently high (i.e., order of 10^ cm"3) and the electron 
temperature is sufficiently low (i.e*, approximately 1200 °k) so 
that electron-ion recombination* is the dominant electron loss 
mechanism . 

The voltage-current characteristic in the abnormal glow 
region is shown in Figure 2 for various pressures. At low cur- 
rents, there is a negative resistance phenomenon which is not 
displayed in these curves. These characteristics, as well as the 
plasma parameters, are dependent upon the helium purity and the 
condition of the cathode surface. Since a large volume exists 
with a large aluminum cathode and significant cathode and anode 
dissipation, the gas purity and cathode surface cleanliness are 
dependent upon the discharge current and the helium flow rate. 

No attempt was made to obtain an ultraclean cathode since, due 
to its size, this is nearly impossible. This discharge is 
intended as a means for generating steady-state, well-behaved 
plasmas of large density and physical size* Therefore, the 


*This is most likely collisional-radiatlve recombination 
(ref. 3) with, perhaps, some disassociative recombination at the 
higher pressures of operation. 
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Figure 2 Voltage-current characteristics of the V-groove cathode discharge 
at various pressures. 








reproducibility of the parameters is important. This repro- 
ducibility was obtained by operating the discharge with the same 
helium flow rate at each pressure, and only recording data at 
each current after the discharge had stabilized. With this pro- 
cedure, excellent data reproducibility was obtained. 

From Figure 2, it is noted that the operating voltage is 
substantially lower than that of the brush cathode discharge 
developed by Persson (ref. 2). This is due to the higher second- 
ary emission coefficients of aluminium compared to these co- 
efficients for the materials used by Persson. Thus, not only is 
the V-groove cathode easier to fabricate,* but it also is subject 
to far less cathode dissipation for the same discharge current. 
Furthermore, the lower beam voltage results in a larger source 
function for electron-ion pair generation near the cathode, 
because the ionization mean free path is smaller at these lower 
beam voltages. Typically, the voltage results in beam energies 
yielding electron-ion pair source functions near the maximum 
value in helium. Also, the lower operating voltage tends to 
lessen the arcing which sometimes occurs near the glass container 
at the cathode. However, the disadvantage of the lower operating 
voltage is that the primary beam- reaching distance is smaller 
than that of the conventional brush cathode discharge, and since 
electron-ion recombination is the dominant electron loss mechanism 
this results in a substantial decay in the secondary electron num- 
ber density in the axial direction away from the cathode. This 
effect occurs for all discharge currents up to 900 mA at pressures 
above 0.8 torr and at sufficiently low voltage and current for 
pressures below this value. The beam-reaching distance can be 
calculated using (ref. 2) 


L — _JL_ I dv (1) 

PV, J (V) 1 } 

1 0 

where $ (V) is the inverse ionization mean free path, S is a mea- 
sure of the ionization efficiency, p is the pressure in torr, 
is the ionization potential of helium, and V is the cathode 
fall voltage which is assumed to be equal to the discharge voltage 
Using 3 = 3.0 and assuming <fi{V) = C/V, where C - 430 in helium, 
the beam-reaching distances for various currents and pressures 
are shown in Table I. Only for beam-reaching distances less than 
15 cm is a substantial effect of beam decay on the axial varia- 
tion of secondary electron number density expected. Results 


*This is of particular importance for this experiment since 
the cathode is so large. 
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which are consistent with this argument will be presented later. 
These results assume: 

(1) g is constant with electron energy; 

(2) The ionization mean free path varies inversely with 
electron energy; 

(3) The beam is monoenergetic at the cathode. 

Assumption (1) is obviously incorrect but, since 3 is not 
known accurately at any energy except that it lies between 1 
and 10, the use of 3 » 3.0 for all energies appears to be rea- 
sonable for the qualitative argument which is presented. Assump- 
tions (2) and (3) are also questionable at the low beam voltages 
used in this experiment. Hence, the beam- reaching distances 
shown in Table I are, at best, crude approximations but they will 
be used to support the hypothesis that the V-groove cathode re- 
sults in a beam-generated plasma. 

Finally, the visual appearance of the discharge is described. 
At applied voltages comparable to the normal cathode fall volt- 
age, a weak blue-green glow appears in the chamber. As the volt- 
age is increased above this value, the glow quickly changes to a 
deep pink color in the center of the discharge reqion and is 
surrounded by a blue-green boundary near the electrodes and the 
glass wall. The thickness of this boundary decreases with in- 
creasing pressure. The pink color is indicative of recombination 
in helium, and the blue-green boundary defines the region where 
the primary electron loss mechanism changes from volume recombina- 
tion to ambipolar diffusion to the walls. For small operating 
voltages at the lower pressures (i.e., 0.4 to 0.8 torr) and for 
all operating voltages at higher pressures, the glow is most 
intense near the cathode. At high operating voltages and cur- 
rents and at pressures less than 0.8 torr, the discharge is visu- 
ally symmetrical relative to the center of the chamber. 


MEASUREMENTS OF THE PLASMA PARAMETERS 


A movable Langmuir probe was used to obtain electron density 
and temperature as a function of position, discharge current, and 
gas pressure. For the range of these parameters over which the 
discharge was operated, the electron temperature is constant. A 
typical set of Langmuir probe current-voltage curves at a fixed 
pressure (0.6 torr) and a fixed position (near the center of the 
discharge chamber) for discharge currents ranging from 0.1 to 0.9 
amp is shown in Figure 3. Since the saturation electron current 
can be determined quite accurately from these data and the electron 
temperature is constant, the reduction of these curves to semi log 


7 



Figure 3.- Typical Langmuir probe current-voltage curves for a 
probe located at the center of the discharge chamber 
for p = 0.6 torr and for discharge currents from 
0 . 1 to 0 . 9 amp 

curves was not done in the bulk of the data reduction. Figures 4 
through 10 show electron density profiles in the axial direction 
for various pressures and discharge currents. The independent 
variable is normalized to 3 inches and the electron density is 
normalized to its peak value. Figures 11 through 17 show the 
peak electron density as a function of discharge current for var- 
ious pressures. The radial electron density profile was also 
obtained but is not of much concern to the results of this invest 
igation. Previous investigations (refs. 2, 4) have shown that 
Langmuir probe data are questionable in a beam-generated plasma. 
Wo such difficulty was encountered with the V-groove cathode ex- 
cept when the probe was placed quite near the cathode (i.e., with 
in the blue-green boundary near the electrode) . It is believed 
that this is due to the lower operating voltage of the discharge 
which gives a less energetic beam. 
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NORMALIZED ELECTRON DENSITY 


P = 0.3 TORH 



NORMALIZED AXIAL DISTANCE 


Electron density profiles for p = 0.3 torr and 
discharge currents ranging from 0.1 to 0.9 amp; 
comparison of theoretical curves and experimental 
points . 


Figure 4 
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Electron density profiles for p - 0,4 torr and 
discharge currents ranging from 0,1 to 0,9 amp; 
comparison of theoretical curves and experimental 
poin 
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Electron density profiles for p = 0.6 torr and 
discharge currents ranging from 0.1 to 0.9 amp; 
comparison of theoretical curves and experimental 
points . 
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Figure 9.- Electron density profiles for p = 0.8 torr and 
discharge currents ranging from 0.1 to 0.9 amp; 
comparison of theoretical curves and experimental 
points . 
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Figure 10.- 


Electron density profiles for p = 0.9 torr and 
discharge currents ranging from 0.1 to 0.9 amp; 
comparison of theoretical curves and experimental 
points . 
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Figure 12.- Comparison of the peak electron densities measured 
by the Langmuir probe and the electron densities 
obtained from the ref lectometer as a function of 
discharge current for p = 0.4 torr. 
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Figure 13. - Comparison of the peak electron densities measured 
by the Langmuir probe and the electron densities 
obtained from the ref lectometer as a function of 
discharge current for p = 0.5 torr. 
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Figure 14.- Comparison of the peak electron densities measured 
by the Langmuir probe and the electron densities 
obtained from the ref lactometer as a function of 
discharge current for p = 0.6 torr . 
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ELECTRON DENSITY 
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Figure 15.- Comparison of the peak electron densities measured 
by the Langmuir probe and the electron densities 
obtained from the ref lectometer as a function of 
discharge current for p = 0.7 torr. 
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Figure 16.- Comparison of the peak electron densities measured 
by the Langmuir probe and the electron densities 
obtained from the ref lactometer as a function of 
discharge current for p — 0.8 torr . 
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Figure 17.- Comparison of the peak electron densities measured 
by the Langmuir probe and the electron densities 
obtained from the ref lectometer as a function of 
discharge current for p = 0.9 torr. 
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An X-band microwave ref lectometer was connected to an open- 
ended waveguide port positioned in the center of the anode. The 
reflection coefficient of the open-ended guide plus plasma layer 
was monitored as a function of discharge current for fixed fre- 
quency and as a function of frequency for fixed discharge current. 
A typical curve of reflection coefficient as a function of dis- 
charge current is shown in Figure 18. In either case, the 
transition from low to high reflection is well defined and de- 
termines quite accurately the condition Up = oi, where uip is the 
plasma frequency and w is the wave frequency of the reflecto- 
meter. Except for operating conditions where there is a sizable 
beam decay, the plasma frequency is described by the peak 
electron density. The electron densities obtained by this 
method are compared with the peak densities measured by the 
probe in Figures 11 through 17 . The agreement is consistent 
with the accuracies of the two measurement techniques. 

There is no accurate method for measuring the collision fre- 
quency of a large, dense plasma. However, from knowledge of the 
electron temperature, the gas type and pressure, and the ion 
density, it is possible to calculate the collision frequency. 

The electron-neutral collision frequency is determined from pub- 
lished data on the collision cross-section for momentum transfer 
in helium (ref. 5) using 


v 


en 


lOOp P 



( 2 ) 


where p is the pressure in torr, P c is the collision probability 
in inverse cm at 1 torr, k is Boltzmann's constant, T e is the 
electron temperature, and m e is the electron mass. For the mea- 
sured parameters of the discharge, v en ranges from about 1.5x10 s 
to 4.5xl0 8 for pressures between 0 . 3 ~and 0.9 torr. At the rela- 
tively low electron temperature and high electron densities gen- 
erated by the V-groove cathode discharge, electron-ion collisions 
are expected to be important (refs. 2, 6). To verify this, the 
following expression for the electron ion collision frequency is 
used (ref . 5 ) : 


1/2 in <9n D > 
V ei W 9n D 


(3) 


where lu is the plasma frequency and n D is the number of electrons 
in a Debye sphere. For electron densities corresponding to X-band 
cut-off frequencies and the temperature of this experiment, this 
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expression gives values larger than the electron neutral colli- 
sion frequency by a factor of about 2 to 3 depending on the 
pressure. Of course, it is questionable whether expression (3) 
is accurate to this factor. In any event, it is expected that 
\> e ^ is at least comparable to v en . The total collision frequency 
is 


v t “ v en + 


V 


ei 


(4) 


THE ONE- DIMENSIONAL DIFFUSION- RECOMBINATION THEORY 

Qualitative arguments have been presented to show that the 
plasma is beam generated. The electrons in this beam emanate 
from the cathode and are rapidly accelerated in the cathode fall 
region. This region is assumed to be very short compared to the 
anode-cathode separation? that is, the beam electrons leave the 
cathode with energy corresponding to the full cathode fall volt- 
age which is assumed equal to the discharge voltage. The beam 
electrons produce the plasma electrons and ions by impact ioniza- 
tion.* The expression <f» (V) is defined as the inverse ionization 
mean free path at 1 torr for an electron of energy V, so that 
the number of electron- ion pairs generated per second per unit 
volume, hereafter called the source function for electron-ion 
pairs, is, 


S 



(5) 


where J is the beam current density which is assumed equal to the 
discharge current density, p is the discharge pressure in torr, 
and e is the electron charge. The beam energy at any axial point, 
z, is the initial beam energy minus the ionization and excitation 
energies expended in reaching the point in question. A one- 
dimensional planar problem which neglects variations in the 
radial direction will be considered. This is a good approxima- 
tion for determining the axial electron density profile, since 
the cylinder diameter is twice the electrode spacing. 

It is difficult to obtain a good functional representation 
for S{z) in terms of the discharge voltage, current, and position 
z. Instead, consider the source function to be known and 
some idealized cases are treated later. If electrons and ions 


*It can easily be shown that the plasma density is orders of 
magnitude greater than the beam electron density. 
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are lost by re combi nation and ambipolar diffusion, then, in the 
steady state: 


D — S- - an 2 + S(z) = 0, 
dz^ 


-d <, z <. d 


(6) 


where D is the ambipolar diffusion coefficient, a is the recombi- 
nation coefficient, and n represents either the electron or ion 
number density. Both D and a are assumed to be independent of 
position and the following boundary conditions are used: 


n (d) = n (-d) = 0 . 


Equation (6) may be written: 

+ 6 2 [s£x) - yf 2 ] = 0 , (7) 

dx 

where 

f = n/n Q , 
x * z/d, 


S (x) = S q s (x) , 



and n Q is the peak electron density and S D = S(o). 

First consider the special case s (x) = 1. This is a good 
approximation to the experimental situation when the beam- 
reaching distance is large compared to the distance 2d. Thus, 
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the source function is constant and the solution is symmetrical 
about x * 0 , where the electron and ion densities are maximum. 

With the substitution y = 6x, the following expression is obtained. 


+ 1 - yf 2 = 0, 0 < y < 6 , 

dy Z 


( 8 ) 


subject to the boundary conditions 


d. f 

f(0) = 1, If 


= 0, f [6 ) - 0 


x=0 


Equation (8) has the first integral 


Ki) + f - 1 - K* 3 - 0 - « 


(9) 


Using the substitution = 1 - f, the solution of Eq. (9) is 
written as 


-Vlf -p=^= 


(10) 


3C 2 - C 4 


which can be expressed in the form of an elliptic function of the 
first kind (ref. 7) F(0\8). Thus 


y = 


VI 

n f(6 \ 8) 


( 11 ) 


1 + 37 (1 - Y) J 


where 


tan 8 = 


V 


3? (1 - Y) 


V 


( 12 ) 


1 + ^ (1 - 7 ) 


- 1 
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V 1 ”” + 3? (1 * y) g 4 


_1 + Vw (1 " Y, ][V 1 + w 11 ~ 


y) - 1 + | g 2 


5 


. ( 13 ) 


The quantity 6 is obtained from 


6 = 


V? 




1 + 37 (1 


- Y) 


-=TX7% F(e l N 8) 

J 


(14) 


where 


sin 0i = 


V 1 + 5 u-u 


i + Jyj- (l- y) l + 3^ 


(1- y) ||1 + 3t/l + 37 (1- Y) 




. (15) 


Figure 19 shows normalized electron density as a function of 
normalized axial distance for various values of y. The value 
y = o is the zero recombination limit. The limit of zero diffu- 
sion where all electrons are lost by recombination is y = 1, 
and <5 = °°, where f = 1, for 0 <. x < 1. From Figure 19, it is 
noted that even when y is quite near 1 so that nearly all the 
electron loss is by recombination, the profile is still greatly 
influenced by diffusion. 

When the electron-ion pair source function is not constant, 
the functional behavior of S(z) must be known to obtain solutions 
of Eq. (6). However, there is probably no one functional form of 
S (z) which is valid for all discharge voltages and currents. In 
order to display in a qualitative fashion that beam attenuation 
can account for the asymmetrical profiles which occur at suffi- 
ciently low voltage or sufficiently high pressure, the second- 
order, non-linear differential Eq. (6) will be solved with 

S (z) = S 0 e" z/L . (16) 
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NORMALIZED AXIAL DISTANCE 


Figure 19.- Theoretical electron density profiles showing 
normalized electron density as a function of 
normalized axial distance for various values of 7 . 


29 


Hence Eq. (7) becomes 


sif + 6 2r e -(d/L) X . yf 2] . 0 f 

dx^ 


(17) 


and solutions of this equation must be obtained numerically. To 
facilitate the numerical procedure, substitute 


y = & 1 (x + x x ) 


where x^ is defined such that 


Hence 


dn 

dx 


= 0, n(-x^) = 


n. 


X=~X- 


dff 

dz 2 


+ e 


-Ry _ 


Yj_f = 0 , 


(18) 


where 


- (d/L)x, 


Y 1 = Ye 


2 2 (d//L)x l 

6 jL = 6 e 


R = 


5 1 L ' 


and the boundary conditions are 


f [5]_ (1 + xj)] = f[6 1 (x 1 - 1)] = 0 
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For given R and y lf Eq. (18) can be integrated numerically 
for both positive and negative y using the initial conditions 
f{0) ® 1, f'(0) = 0, and the boundary conditions yield the con- 
stants 6-^ and X]_. These constants are then used to yield solu- 
tions in terms of the normalized variable x. 

The procedure outlined above for treating the case of a 
non-constant source function simplifies the numerical computa- 
tions for finite diffusion. However, the attenuation constant R, 
which is specified in this numerical procedure, only gives the 
actual constant L, once Sg is determined (i.e., L = d/R6x) . Also, 
the recombination limit is not obtained for finite R. That is, 
the recombination limit for finite L corresponds to 


<5 *- oo 



Y — ► 1 


and hence 


R —*► 0 



such that R<5]_ is finite. In this limit, yg — ►exp (~R6]_) < 1. 
However, the numerical procedure is greatly simplified using fixed 
R and experimentally the full recombination limit is never at- 
tained. A best-fit procedure was employed to compare experiment 
and theory using a sizable number of theoretical curves for wide 
ranges of and R, including R = 0. These comparisons are 
shown in Figures 4 through 10. The notations y and 6 are used in 
these figures for Yy and 6p. In view of the experimental error 
of the Langmuir probe, the inadequacy of the exponential form 
for approximating the source function and the finite amount of 
computer data used for matching, these comparisons are made only 
to support the basic contentions that the plasma is beam-genera- 
ted and that recombination plays a dominant role in determining 
the electron densities. Values of L computed from the parameters 
R, y, and 6 for best matching are comparable to the beam- reaching 
distances computed in the second section (except, of course, for 
cases where R = 0 is assumed) , but there is no exact correspond- 
ence in magnitude, or even in relative values, with changes in 
either pressure or discharge voltage. 
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Finally, the recombination coefficient will be estimated 
from the measured discharge parameters. For cases where R •» 0, 
the values of y which give the better match between theory and 
experiment are all close to 1. Thus, for a small volume near 
the center of the discharge, virtually all electron loss is by 
recombination and 


n 


o 


J eM m 

yf ea 


/ 


or 


a 


JpS (V) 

2 * 


en o 


Assuming V to be the discharge voltage and J to be the discharge 
current density, a is calculated. Values of a obtained in this 
manner are shown in Table I. These values are consistent with 
the collisional-radiative recombination theory (ref. 3). This 
theory uses the hydrogen atom, but results for helium should not 
be significantly different. In any event, it was concluded that 
the discharge is predominantly recombination-controlled for opera- 
ting conditions such that R =» 0 . No inference should be drawn 
from the fact that a appears to increase with either decreasing 
discharge voltage or increasing pressure. This tendency can more 
likely be attributed to the assumptions concerning the equality 
of beam and discharge parameters (i.e., voltage and current). 

As the beam voltage is decreased or the pressure is increased, 
the contribution of ion flow to the total discharge current is 
expected to be increasingly important. Also, under these condi- 
tions the beam has an increasing spread of electron energies. 


CONCLUSIONS 

The results of an experimental investigation of the plasma 
parameters generated by a large V-groove cathode have been pre- 
sented, The plasma can be described by the following parameters 


n =» 2xl0 12 cnT 3 


d a 15 cm 


v ~ IQ 8 - 10 9 sec 1 


T = 1200°K 
e 
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Furthermore, the plasma density is recombination-controlled, and 
electron-ion collisions play an important role in determining the 
total collision frequency for momentum transfer. 

The results described in this technical note have been con- 
fined to discharge currents less than 1.0 A and an electrode 
spacing of 6 inches. The discharge has recently been operated 
with currents up to 5 A and an electrode spacing of 12 inches. 
Complete diagnostics are not yet available, but microwave inter- 
ferometer measurements and ref lectometer measurements indicate 
electron densities greater than 4 x 10^2 e/cm^. 


Electronics Research Center 

National Aeronautics and Space Administration 
Cambridge, Massachusetts, March 1969 
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